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Abstract 

A new type of basis functions is proposed to describe a two-electron continuum which arises 
as a final state in electron-impact ionization and double photoionization of atomic systems. We 
name these functions, which are calculated in terms of the recently introduced Quasi Sturmian 
functions, Convoluted Quasi Sturmian functions (CQS). By construction, the CQS functions look 
asymptotically like a six-dimensional spherical wave. The driven equation describing an (e, 3e) 
process on helium in the framework of the Temkin-Poet model has been solved numerically using 
expansions on the basis CQS functions. The convergence behavior of the solution has been examined 
as the size of the basis has been increased. The calculations show that the convergence rate is 
significantly improved by introducing a phase factor corresponding the electron-electron interaction 
into the basis functions. Such a modification of the boundary conditions leads to appreciable change 
in the magnitude of the solution. 
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I. INTRODUCTION 


The Coulomb three-body scattering problem is one of the most fundamental outstanding 
problems in theoretical atomic and molecular physics. The primary difficulty in description 
of three charged particles in the continuum is imposing appropriate asymptotic behaviors of 
the wave function. 

Several ab initio methods are developed for constructing solutions to the three-body 


scattering problem (see the review [1]). The exterior complex scaling (ECS) method (see 2| 
and references therein) allows the problem to be solved without explicit use of the asymptotic 
boundary conditions. Specifically, ECS recasts the original problem into a boundary problem 
with zero boundaryconditions. (For an extension of ECS to the case of long-range Coulombic 
interactions see [ 3 , 4|.) Some of the other methods use a product of two fixed charge Coulomb 
waves to approximate the asymptotic three-body continuum state. The convergent close 
couplin g (C CC) @-.7] and the Coulomb-Sturmian separable expansion B a and the J- 
matrix 10j, [ll| methods treat the problem in the Laguerre basis representation. The latter 
two methods transform the original problem to a Lippmarm-Schwinger-type integral equation 
whose kernel seems to be generally non-compact. Alternatively, the Generalized Sturmian 
Functions (GSF) method 12, |l3| converts the problem into an inhomogeneous Schrodinger 
equation with a square integrable driven term. One-particle generalized Sturmian functions 
with an appropriate asymptotic behavior are obtained (numerically) as eigensolutions of a 
Sturm-Liouville problem. The GSF method driven equation is solved by an expansion into a 
basis set of two-particle functions, which are products of two generalized Sturmian functions 
that both satisfy outgoing-wave boundary conditions fu3 |. 

In the present paper in order to describe a Coulomb three-body system continuum we 
propose a set of two-particle functions, which are calculated by using recently introduced 
so called Quasi Sturmian (QS) functions jl5j. The latter satisfy a two-body inhomogeneous 
Schrodinger equation with a Coulomb potential and an outgoing-wave boundary condi¬ 
tion. Specifically, the two-particle basis functions are obtained, by analogy with the Green’s 
function of two non-interacting hydrogenic atomic systems, as a convolution integral of two 
one-particle QS functions. The QS functions have the merit that they are expressed in closed 
form, which allows us to find an appropriate integration path that is useful for numerical 
calculations of such an integral representation. We name these basis functions Convoluted 
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Quasi Sturmian (CQS) functions. Note that by construction, the CQS function (unlike a 
simple product of two one-particle ones) looks asymptotically (as the hyperradius p —>■ oo) 
like a six-dimensional outgoing spherical wave. 

We apply these CQS functions to the solution of a problem of double ionization of He in 
the framework of the Temkin-Poet model. We solve the driven equation describing an (e, 3e) 
process jl4] by an expansion into the basis set of CQS functions and explore the convergence 
properties of the expansion. Note that the CQS functions asymptotic behavior in the so 
called three-body region Oo where all three particles are well separated is not correct since 
it misses out the phase factor, corresponding to the Coulomb interelectronic interaction. 
Therefore, the expansion method effectiveness is open to question. In order to improve the 
convergence rate, we equip the basis functions with the phase factor corresponding to the 
potential —. 

The paper is arranged as follows. In Sec. II we present the three-body driven equations 
jldl] and [l6| whose solutions possess all the information for the (e, 3e) process on helium and 
that for the one-photon ionization, respectively. Here we suggest the CQS functions which 
form a basis set used for solving these equations. The CQS functions are expanded in a series 
of products of the single-particle Laguerre basis functions. A useful integral representation 
is also introduced for the CQS functions. The asymptotic behavior of the basis functions in 
the region Q 0 is deduced from their integral representation. In this section we also propose 
a modification of the basis functions which allows to take into account the e — e interaction. 
The use of both original and modified versions of the basis functions in solving the s-wave 

I f 

driven equation |14j is considered in Sec III. Finally, Sec. IV provides a summary. Atomic 
units are assumed throughout. 


II. QUASI STURMIAN BASIS FUNCTIONS 
A. Driven equations 

Electron-impact ionization and double photoionization of atomic systems can be cast 
as an inhomogeneous three-body Schrodinger equation with a square integrable right hand 
side. For example, in the approach {jjJ to the (e, 3e) process on helium, the four-body 
Schrodinger equation is reduced to the following driven equation for the three-body system 
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(e ,e , He ++ ) = (1,2,3): 


E-H 


$ (+) (ri,r 2 ) = W / /i (r 1 ,r 2 )<l) (0) (ri,r 2 ). 


( 1 ) 


k 2 k 2 . 

E = -j- + -j- is the energy of the two ejected electrons. The three-body helium Hamiltonian 
is given by 


h = h 1 + h 2 + —, 

r 12 

1 2 

Hj = —-A r ,-, J = 1,2. 


( 2 ) 

( 3 ) 


$(°)( r i, r2 ) represents the ground state of the helium atom. The perturbation operator Wfi 
is written as 

^ 1 A'TT 

( 4 ) 


W7i(ri, r 2 ) = ^^(-2 + e* qri + e^ 2 ), 


(27t) 3 q 2 

where q = kj — kj is the transferred momentum, k ; and k f are the momenta of the incident 
and scattered electrons. 

In turn, the one-photon ionization problem also takes the form of the driven equation 
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E-H 


$ (+) (ri,r 2 ) = -So ■ z3 G $ (0) (n,r 2 ), 


( 5 ) 


where So is the amplitude of the electric-held vector and Dq is the dipole operator. 


B. Convoluted Quasi Sturmians 


Our method of solving the driven equations (JT]) and (J5|) is to expand the solution 


N—l 


4> l+) (ri.r 2 ) = E E c S 2) I«i<W2;£M> 


Q 


L,C,\ n i,ri2=0 


on the basis 


KW 2 ;£M) Q = 


nr 2 


y £ LM ih,h)= (^i m i^rn 2 \LM)Y hmi (r l )Y hm2 (r 2 ). 

17111712= M 

Each function Qnfni + ^ is assumed to satisfy the radial equation 

Vn\( r l)Vn 2 ( r ^ 


E-h['-hf g« 2 2 W(E;r 1 ,r 2 ) = 


7l\ V ' 77.2 ' 

r\r 2 


( 6 ) 

( 7 ) 

( 8 ) 


( 9 ) 
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where 


( 10 ) 


ll = - 


1 d 2 1 1(1 + 1 ) 2 

2 <9r 2 2 r 2 r ’ 


are the Laguerre basis functions (b is a real scale parameter) 


^n(' r ) = i( n + T-Wi] i (2hr) £+1 e br L 2 n i+1 (2br), 


( 11 ) 


which are orthogonal with the weight 


l. 


drtfjr)-tf m {r) = S nm- 


( 12 ) 


In order to obtain the Qn-fni^ with the outgoing-wave boundary condition we use the 
Green’s function G''' p/2 ' + > (E) which can be expressed in the form of the convolution integral 

13, U 

QhH+)f E ) = — [ d£G h{+ \V2£)G^ + \J2(E - £)), (13) 

2tti J 
c 

where the path of integration C in the complex energy plane £ runs slightly above the branch 
cut and bound-states poles of (j h+) (see Fig. [[]). Applying the Green’s function operator 
(which is the inverse of the operator in the left-hand side of ()9j) ) onto both sides of equation 
dnj), we find that 


QZ% +) (E-,r u r 2 ) = — / d£ Q*< + >(t/2 £■ n)Q« + > (yj2(E - £); r 2 ), (14) 


where the one-particle quasi Sturmian finictioii Qn is defined by 15] 


Q^(fc,r)= / dr'G f(±) (l;r,r')—tZ'n(r'). 


(15) 


We name the basis functions (1131) Convoluted Quasi Sturmian (CQS). 


C. Laguerre basis expansion 

CQS functions can be expanded in terms of the Laguerre basis functions (fill) as 

<&™ <+) (E;n,r a )= £ (16) 

mi,7712=0 
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The coefficients G'mimliin 2 (E) are the matrix elements of the Green’s function (1131) over the 
functions 

V’niOr)^^) 


and can be calculated using the convolution |9 
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/ MC%£(V2£)G‘*»<.s/2(E-£)) 


(17) 


(18) 


Cl 


of two one-particle Green’s function G matrix elements. The latter is expressed in terms 


of two independent J-matrix solutions 
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20 ]: 


S„,(k) = 1 [(n+ l) (2t+1 ,] 1/2 (2 sin g) w e-W 2 or* |r( gh)f 1 

x(—t u) n 2 T 1 (—n, ■£ + 1 + i/3; 2£ + 2; 1 — cu 2 ), 

C&’M = -i/«!(« + 2f+l)^^ 

x ihJh; r ( —)( j i-K'i -k (-1 + iP,n + l;n + 1 + 2 + iftu 2 ), 


where 


w = e* = 


6 + ik 


sin£ = 


26fc 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


b — ik' b 2 + k 2 ' 

(3 = ^ is the Sommerfeld parameter. Following the method of |l8|, we perform the inte¬ 
gration in (gB) along contour C\ (see Fig. [T]) which is obtained by rotating the contour C by 
angle —7r < <p < 0 about the point This allows us to avoid the singularities of Gm[n\- 
Figs. |3] and [4] show the behavior of Qq^ +) (TTGI) (the dashed lines) along the = r 2 = p/V2 
diagonal. For these calculations, we choose E = 0.735 a.u. and 25 for the upper limit of the 
sum. The scale parameter b was set to 1.6875 and the rotation angle <p was — 


D. Asymptotic Behavior 


The asymptotic form of the QS function Qn ± ' ) can be written as 15 1 


where 


Qi {± \k,r) - 1 S nl {k ) e ±i(kr-fnn(2kr)-f +*/(*)) ? 

i—>oo /c 


0 i<Ti.(k) _ r(i? +1 + i(3) 


\r(£ + i + i/ 


( 23 ) 
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The asymptotic behavior of the CQS function (fT4l) for r\ —> oo and r 2 —> oo simulta¬ 
neously (in the constant ratio tan(a) = r 2 /ri, where a is the hyperangle) is obtained by 
replacing Qni +} and Qn} + ^ by their asymptotic approximation (1231) and making use of the 
stationary phase method to evaluate the resulting integral. The stationary point £ 0 which 
satisfies the equation (see, e. g., 2l|) 


— (v^rq + VE*Er^ = 0. 
is £ 0 = cos 2 {a)E. Therefore, we finally obtain 


(24) 


Q 


nin 2 


(E]n,r 2 ) 


e*S nitl ( Pi)S n2 e 2 (p 2 )4 


p—too •“ v “ ‘ y/~P 

x exp <; i y/2Ep - A ln(2pin) - A ln(2p 2 r 2 ) + a tl (pi) + a ta (jp 2 ) - 

where p = x/rf+rf is the hyper-radius, pi = cos (a)V2E, p 2 = sin(a)\/2 ~E, A ,2 = 


(25) 


~2 
Pi,2 ' 


E. Integral representation 

Rather than apply the expansion (ORl) to calculate the CQS functions at large distances, it 
might be more convenient to employ the contour integral (1T41) whose integrand is expressed 
in terms of the integral 0 


Qn {± \k,r) = - [(n + 1 ) 21 + 1 ] 2 (2 br) 


£+1 g—br 


wm f (i - *)““(! - 

°. ( 26 ) 

x (1 — z — u> ±1 z) n exp (z [b ± ik] r) L^ £+1 ^ ~~~rr ^ 2 brj . 

Note that a part of the rotated straight-line contour C\ indicated by a dashed line in Fig. [T] 
lies in the unphysical energy sheet (—27r < arg(£) < 0). In turn, Qn + \y/2S,r) diverges 
exponentially for large Im(£) (see, e.g., ([23]) ) in the lower half-plane. In order to ensure 
convergence of the integral (TT4l) we deform the contour C\ in such a way that the resulting 
path C 2 , shown in Fig. [2j asymptotically approaches the real axis. Specifically, the energy S 
on the contour C 2 is parametrized in the form 

.(!-*) 


S — t + iD- 


(27) 


l + ( 2 ’ 

where D is a positive constant and t runs from 00 to — 00 . Generally there are points on 
the contour C 2 at which Re(f? + i/3) < —1 and thus the integrand in (1261) is singular at the 
endpoint z — 1. To avoid this, we apply the following procedure. Let m be the minimum 
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positive integer number such that — m < Re(£ + i/3). Then integrating (1261) by part m — 1 
times and assuming the integrated terms vanish at the limit z = 1 we obtain the integral 
which can be evaluated numerically. As a test of this analytic continuation of Qn + \k, r) 
into the complex fc-plane, the CQS function Qqq^ has been calculated using the integral 
dH along the contour C 2 . We choose D = 0.85 for which m = 3. The results for real and 
imaginary parts are plotted in Fig. [3] and Fig. [4] by the solid lines. We also plot the CQS 
function Qqq calculated with D — 15 for which —1 < Re(i/3), i.e., rn — 1 (the dashed 
lines). Agreement between these representations of CQS functions and (TT6l) argues for the 
suggested analytic continuation. 

Note that the asymptotic behavior (1751) of the two-particle CQS functions depends upon 
the indices ri\ and n 2 ■ It follows from (12U1) that this dependence can be eliminated by 
dividing (TT4j) by (p 1 )R r ( 2 2 (p 2 ), where 

Bn(k) = [(n + 1) 2 ^+i] 1//2 (—o;) n 2 Fi (— n, £ + 1 + i/3; 21 + 2; 1 — uj 2 ) . (28) 

We present in Figs. [5] and [6] a few s-wave CQS functions /B^ l (p 1 )B^ 2 (p 2 ) for a — |. 

These functions asymptotic behavior at large distances is shown in Figs. [7] and [8j For 
comparison, we also show the asymptotic approximation (|25l) for Qqq^. 

Note that it follows from (1261) that on the left part of the contour C 2 where k ~ i\k\ and 
\k\ —y 00 the function Qn for large r behaves like e~ br rather than e lkr . Thus, the greater 
is the scale parameter b, the faster the CQS function (fl4l) reaches its asymptotic form (|25]h 


F. The solution asymptotic form 


We try to solve the equation (pQ) by an expansion into the basis set of CQS functions (1T|) 


whose asymptotic behavior in the region fl 0 is n °f correct since it misses out at 


eas 


t the 


phase factor, corresponding to the Coulomb interelectronic interaction (see, e. g., 21, 22]): 

P 1 


W 3 ( r i, r 2 ) = 


V2Er V2 


In ( 2 V 2 ~Ep 


(29) 
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Inserting (1251) into © we find the formal result for the the asymptotic form of the solution 
of © at large distances: 


$ (+) (ri,r 2 ) 


ivr exp©\/2i?p— «i ln(2piri )—02 ln( 2 p 2 C 2 )] } 

e 4 ^72 


Esin(2(f)) Y 7r 

x E y e LM (ri, r 2 ) exp <j i a h (p i) + a h (p 2 ) - ^ (£l 2 +fa) 


W 2 n 


(30) 


iV-l 


X X) C'n 1 ( n 1 2 £2) ^n 1 £ 1 (Pl)-?n2r 2 (P2) 
n 1 ,n 2 =0 


Thus the feasibility of using the CQS functions © for solving © is questionable. To 
improve the asymptotic properties of the CQS function (and thereby solve the problem of 
slow convergence of the expansion ©), it may be useful to modify the CQS function by 
multiplying it by 


e *W V2 (ri,r 2 ) 


where 


Whh (r 1 ,r 2 ) 


P 


^ v ^V (l& (r 1 ,r 2 )l, 1 (2v / 2Ip), 


1 


V^ 2 (ri,r 2 )= / dfidf 2 Km (h, r 2 )] *— r 2 ). 

J 02 


Hence the modified function takes the form 


Q'yjy (S; r!, r 2 ) = {E . ri _ r2 ). 


(31) 

(32) 

(33) 

(34) 


Note, however, that such phase factors can not take into account the off-diagonal elements 

(ri,r 2 ) = f dr l( lr 2 [y 1 ^ (fi, r 2 )] * — Km ( f i, r 2 ). (35) 

J r 12 

Thus, it will probably be more convenient to choose the basis functions of the form 

e iW 3 (r u r 2 ) \ ni l in2 £ 2] LM) q . (36) 

It is beyond the scope of this paper to discuss in detail the modified basis functions and 
their applications. In this report we restrict ourselves to a simple s-wave case. 


III. SOLVING OF THE DRIVEN EQUATION 

In the Temkin-Poet model the equation © reduces to fhlj ] 

Id 2 Id 2 2 2 1 “ 

^ 2 dr 2 2 dr 2 rq r 2 ry 


x(n,r 2 ) = ^(ri,r 2 ), 


(37) 
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where x( r i?^ 2 ) = r i r 2<^~ t4 ( r n r 2 ) ; whereas the right hand side is given by 

J r (ri, r 2 ) = -tttti’T t 2 _ ioO r i) - jo(.qr 2 )] rir 2 — e ~ z ^ ri+r2) . (38) 

{/ 7rj J 7T 

We set E = 0.735, q = 0.24 and Z e = 2 — 5/16. 

To solve the equation we first consider the expansion 


N -1 

x(n,r 2 )= C m,n 2 <3it?n 2 ( ri ’ r2 )- ( 39 ) 

ni,n 2 = 0 

Our discussion is limited to s waves, so that we omit the angular-momentum labels £\ and 
1 2 . We choose Z e for the scale parameter b of the basis. Inserting (TTT71) into (1571) gives 


N-l 

£ 

ni,n 2=0 


^( 0 )^( 0 ) _ 2_ g (+) 
7*17*2 r 


711,712 


(n,r 2 ) 


= ^(ri,r 2 ). 


(40) 


Then, multiplying Eq. (1401) by t/ ; mi(^i)'0m 2 ( r 2) and integrating, in view of the orthogonality 
condition (Tl2|) . gives 


N-l 


^ ^ [^mi,ni^m2,n2 T £mi,m 2 ;rai,n 2 ] f^ni,n 2 77. mijm2 ; r?7 , i, 777 , 2 £ -/V 1, 

rai,ri2=0 


where 


77 r 



dr i dr 2 ^ mi (r 1 )^ m2 Cr 2 )J r (ri,r 2 ). 


(41) 


(42) 


o o 


Further, we approximate the matrix elements 

C 


■'mi, 7712 ;Hl, 712 



0 0 


dridr 2 V’m 1 (7 - 1 )V’m 2 (7 - 2) ( 43 ) 


by using the expansion of CQS function into the Laguerre basis (TT6l) and taking into account 
the basis completeness: 

JV-l 


C 


7711 ,7712 5 711 ,712 


E t/ , , <W+) 

^ mi ,m 2 jri! ,n 2 '- T n ' 1 ,n(, ;m ,n 2 ’ 


(44) 


71^,712=0 


where V miim2 . i7l f n r are the matrix elements of the e — e interaction: 


^TTll ,7712 571^ ,712 



dridr 2 tp mi (ri)V> TO (r 2 )— Vv, MV’njM- 


(45) 


0 0 
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Our aim is to study the convergence properties of the expansion (1391) (in conjunction with 
the approximation (jH|l ) as N is increased. The real and imaginary parts of the solution 
y(ri,r 2 )p 1 / 2 2/sin2a = 0(+)p 5 / 2 along the r 1 = r 2 = p/V 2 diagonal are shown in Figs. 9 and 
10. Figures 11 and 12 show the results for the asymptotic approximation (1HU1) (multiplied 
by p 5 / 2 ) to the solution. From the plots we see that applying the expansion (fTJll yields a 
solution with divergent phase as a function of N, whereas the magnitude 


A = lim \(j)^ + \ri,r 2 )p 5 ^ 2 \ 


p —^OO 


(46) 


seems 


to converge. Actually, from the asymptotic form (1301) of the solution, it follows that 


A n — 


2(2 E) 3 / 4 [2 1 


E sin 2a V n An 


N -1 


Y C ni, na S ni {pi)S n2 {p2) 

ni,ri2=0 


(47) 


i-4 


In the case a — |, pi — p 2 — \[E we have obtained A 16 = 1.505 x 10 4 , A 21 = 1.507 x 10 
and A 26 = 1.400 x 10" 4 . 

It might be expected that the convergence could be achieved using the modified CQS 
functions: 

Q&U’V n) = e ,W(r ‘’ rd Q<il ,(r„ r 2 ), (48) 

P 1 \n(2V2E(l+p)]. (49) 


>V(ri,r 2 ) = 


V2E( 1 + r >) 


Then substituting the expansion 

X(n,r 2 )= Y 


N—l 


ni, ri 2=0 


into (1371) gives 


N—l 


E 

7*1,712=0 L 


^Pni (t'i)V’ 7*2 ( r 2) 


r 1 r 2 


where the operator U is defined by 


u = ~;t + Htr 


1 / 


For large p one hnds from 

d 


that 


i i / aw i _ i 7 a 2 w , a 2 w 
2 ^ ar 2 y 2 _l_ 1F|~ 

aw_a_ , aw a 

dr\ dr\ ' dr 2 dr 2 


Q { nV m (E; rur 2 )~iV 2 E^QM n ,(E-, r u r 2 ), 


dr i jS 


P 


(50) 


n,r 2 ) C„„„ 2 = r 2 ), (51) 


(52) 


(53) 
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and therefore the action of U on Qrt A ln 2 in the asymptotic region is reduced to multiplication 
by the ‘effective potential’ 


U^L(r l ,r 2 ) = 


_ UQ^n 2 (E-, r i, r 2 ) 


ni,ri2 ' 


Q ( nun 2 (E] ri,r 2 ) 


such that 


U nun 2 ( r b r 2) ~ 

’ p—>oo 


In (2y/2Ep 


y/2 Ep 


(54) 


(55) 


Figure [13] shows plots of the and U§1 on the diagonal r\ = r 2 . For comparison we 
have also plotted the e — e potential E 

A matrix equation for the coefficients C ni n2 is obtained by multiplying (15T1) from the left 
by 'ipmi (j' \)ipm 2 (r 2 ) and integrating over both coordinates: 

N—l 


Tl\ ,77-2—0 


^ ^ kmi, ni^m2,ri2 ^ 


m\,rri2\n\,n2 


n \,712 ^mi,m2 5 TYl\. ) TYl2 ^ 1? (56) 


where 


n. 


oo oo 

772 1,7712 J J dr 1 dr 2 ipmi (ri)ip m 2 (r 2 )e 

o o 

iV-l 


£ 


mi ,m 2 ;ni ,n 2 


-iW{n,r 2 ) 7 - 2 ) j 

E tt /'-<(+) 

U mi ,m 2 ,n' 2 Or n , ^ , r 


;ni,n 2 ’ 


n^,722=0 


Um\ ,m2',n' 1 ,n f 2 



dri ir 2 p mi (r 1 )i/j m3 (r 2 )lfl/’„;(r 1 )i/j„. (r 2 ). 


(57) 

(58) 

(59) 


0 0 


Figures IT41 and [T5l show the solution x(r 1; r 2 )p 1//2 2/ sin2o: = 0(+)p 5 / 2 along the diagonal 
ri = r 2 . In turn, hgures [T6] and [TT] present the real and imaginary components of the 
corresponding asymptotic form 

0 ( +)p 5/2 ^ £ 5 ni , n2 ,S ni (p 1 ) < S n2 (p 2 )exp('-f^| s iln(2V2^p)) 

“ m,n 2 =0 v 7 ( 60 ) 

x exp ji V2Ep - ai ln(2piri) - a 2 ln(2 p 2 r 2 ) + cr 0 (pi) + cr 0 (p 2 ) + f } • 

From the plots one can see that the convergence can be achieved by introducing a phase 
factor corresponding to the interelectronic potential l/ri 2 . Furthermore, such a modification 
of the asymptotic form of the basis functions could result in appreciable change in the 
solution magnitude 


A = lim 

p —^OO 


(p i - +) (r 1 ,r 2 )p 5/2 


(61) 
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For comparison, in this case we have obtained A 16 = 7.346 x 10 4 , A 2 1 = 7.396 x 10 4 and 
A 26 = 7.593 x lO" 4 . 

IV. SUMMERY 

Two-particle basis functions — labelled CQS — are proposed. By analogy with the 
Green’s function of two non-interacting hydrogenic atomic systems, they are expressed as a 
convolution integral of two one-particle QS functions. We suggest an analytic continuation of 
the QS functions into the entire complex fc-plane in order to perform the numerical contour 
integration. The asymptotic limit of the CQS basis functions in the region f2 0 is expressed 
in closed form as a six-dimensional outgoing spherical wave. 

We study the application of the expansion into the CQS functions to the solution of the 
inhomogeneous equation describing the double ionization of helium by high-energy electron 
impact in the framework of the Temkin-Poet model. Note that in constructing the basis 
functions we do not take into account the interelectronic interaction, and therefore the 
asymptotic behavior of the basis functions is not correct. Hence, the driven equation whose 
solution is expanded in terms of these basis functions is noncompact (due to the Coulomb 
potential in the left hand side). Thus the applicability of this approach is questionable. 
We show that the problem of slow convergence (or even perhaps lack of convergence) of 
the expansion can be solved by using the modified CQS functions equipped with the phase 
factor corresponding to the potential A--. Moreover, the solutions which satisfy the different 
boundary conditions differ appreciably in magnitude. These results suggest that suitable 
basis functions can be obtained by expanding products of the missing phase factor and CQS 
functions in series of bispherical harmonics. 
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Figure 1: C is the straight-line path of integration of the convolution integral (fl3l) . The bound-state 
poles of (\f 7 2£) are depicted as full circles. The grey line is the unitary branch cut. A part of 

the rotated contour C\ (the dashed line) lies in the region of unphysical energies. 
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Figure 2: The deformed contour C 2 asymptotically approaches the real energy axis. 
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Figure 3: The real part of the s-wave CQS function f° r ^ = 0-^35 and b = 1.6875 along the 

diagonal r\ = r 2 = p/y/ 2, approximated by with the upper limit 25 (dashed line) and obtained 
by integrating (1141) along the contour C 2 (j27|) with D = 0.85 (solid line) and D = 15 (dotted line) 
using the analytic continuation. 
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Figure 4: The same as in Fig. [3] but for the imaginary part. 
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Figure 6: The same as in Fig. [5] but for the imaginary parts. 
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Figure 7: Comparison of the real parts for the first few CQS functions 
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Figure 8: The same as in Fig. [7] but for the imaginary parts. 
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Figure 9: The real components of the solutions ^>(+1 p 5 / 2 for 


different basis sizes along the diagonal 


r\ = r 2 . 
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Figure 10: The same as in Fig. [9] but for the imaginary parts. 
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Figure 11: Real parts of the asymptotic form for the solutions p 5 / 2 for different basis sizes along 
the diagonal r i = r 2 - 
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Figure 12: The same as in Fig. [TT]but for the imaginary parts. 
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Figure 13: Real and imaginary parts of the ‘effective potentials’ Uq q and (1541) . 
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Figure 14: The real components of the solutions for different basis sizes along the diagonal 




Figure 15: The same as in Fig. [14] but for the imaginary parts. 
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Figure 16: Real parts of the asymptotic form for the solutions (f)^ p 5 / 2 (1601) for different basis sizes 
along the diagonal r\ = r?- 
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Figure 17: The same as in Fig. [16] but for the imaginary parts. 
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